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1 Introduction 


Let it be required to solve some system of Diophantine equations. In this case ' alge- 
braic methods can be applied for 

a) constructing the total algebraic solution of the system; 

b) finding the transformations translating an algebraic solution of the system 
from one form into another one; 

c) elucidating the general legitimacies existing between the elements of the alge- 
braic solution; 

d) replacing the total algebraic solution containing L arbitrary selected parame- 
ters by a set of algebraic solutions containing less than L parameters. 

This paper is devoted to further advance of algebraic approach to solving 
systems of Diophantine equations. In particular, in this investigation we 

1) describe the simple way of obtaining a total solution of systems of 
Diophantine linear equations in the integer numbers, and show (see Sect. 2) that this 
way may be considered as some modification of F. Smarandache algorithm 3 from 
his work 2, 

2) demonstrate the effectiveness of the algebraic approach to the elaboration of 
computative algorithms and analytical formulae, which may be used respectively for 
obtaining the required numerical solutions of the discussed systems and for count- 
ing of the total quantity of solutions from a given class of numbers (Sect. 3); 

3) derive analytical formulae available for constructing classical Magic squares 
of both odd and even orders (Sect. 4). 


2 The way of obtaining the total solution of systems of Diophantine 
linear equations in the integer numbers 


Let it be required to solve some system of linear Diophantine equations in the 
integer numbers. It seems to be evident, that there is no complication in solving this 
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problem at present. For instance, one may find in the work? as many as five 
different algorithms to obtain a total solution of this problem, which correctness are 
proved by mathematical methods and illustrated by concrete examples. In particular, 
to illustrate the correctness of an algorithm 3, the system from three following 
equations 


3x, +4x, + 22x, - 8x, =25 
6x, +  46x,-12x,=2 (1) 
4x, +3x, -— x, + 9x, =26 


are solved in the work”. The final solution, obtained for (1) by the mentioned algo- 
nthm, has form 


xX, =—40k, -92k, +27; x, =3k,+3k,+4; x, =—-11k, +8; (2) 
x, = 6k, +12k,-4; x, =3k, -2, 


where k, and k, are any integer numbers. 

Let us clear up a question whether (2) is the total solution of system (1) in the 
integer numbers. To make it we will solve (2) by the algebraic methods with testing 
their correctness on every step of our computations. 


1. As well-known !~3, the total algebraic solution of the system (1) may be found 
by standard algebraic methods (for instance, by Gauss method). In our case it has the 
form 


x, =—(23x, -— 6x, -1)/3 
X, = (x, +2x,+24)/4 (3) 
x, = (-llx,+2)/3 


that coincides with the solution found on the first step of the algorithm 3 of the 
work’. 


2. As well-known from the theory of comparison *‘, the total solution in the inte- 
ger numbers for the last equation of the system (3) has the form x3= {-11(3m + 
mo)+2}/3, where mj, is any integer number; the value of mm is equal 0, +1 or +2 and 
is chosen from the condition that the number (—117 +2)/3 must be integer. Thus, 
we find on second step of our computations that x5= 3m +1 and x3=—11m)-3. 

We note that the solution x5=3k,-2 of (2) may be obtained from our solution 
by change of the variable m,; to &,-1. Thus, both values of xs are identical 
solutions. 

3. Let us get to solving the second equation of the system (3) in the integer num- 
bers. Replacing the value of xs by 3m ,+1 in this equation we obtain that 
X2 = 6 + (6m, + x4+2)/4. Hence it appears (see point 2) that x4= (—2+4]))m)-2+H4h, 
where /; and /, are any integer numbers. 

4. Replacing the value of x4 by (-2+4/,)m),—2+4/, in the first equation of the sys- 
tem (3) we obtain that x) = 2xs—{23(-2+41;)mm; + 921,-47}/3. Hence it appears (see 
points 2 and 3) that J; = 3m+2 and J, = 3m3+1 and, consequently, the total solution 
of the system (1) in the integer numbers has the form 
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xj = —4m, (23m +10)- 92m, — 13; x2 = 3m (mp +1)+ 3773 Bae (4) 
x,=—llm,-3; x, =6m,(2m, +1)+2(6m,+1); x, =3m, +1, 


where m,, m2 and m; are any integer numbers. 

Comparing (4) with (2) we find that 

a) the solution (4) contains greater by one parameter than solution (2); 

b) if m,.=0, m,=k,-1 and m3= k in the solution (4) then the solution (4) coin- 
cides with the solution (2). 

Thus, the solution (4) contains all numerical solutions of the system (1), which 
may be obtained from (2), but a part of numerical solutions, which may be obtained 
by (4), can not obtain from (2) or, in other words, (2) is not the total solution of the 
system (1) in the integer numbers. 

We add that, in general, a partial loss of numerical solutions of systems of linear 
Diophantine equations may have more serious consequences than in the discussed 
case. For instance, as it has been proved in the work°* by using the algebraic 
approach to solving systems of Diophantine equations, 


if Magic squares of 4th order contain in its cells 8 even and 8 odd numbers then 
they can not have structure patterns another than 12 ones, adduced in works Saeed 
for Magic squares, contained integer numbers from I to 16. 


In reality, this statement is incorrect because yet several new structure patterns 
may exist for Magic squares from 8 even and 8 odd numbers *°. 


3 Analysing a system from 8 linear Diophantine equations 


To demonstrate the effectiveness of the algebraic approach to solving some combi- 
natorial and number-theoretic problems, presented in the form of systems of Dio- 
phantine equations, in this section we will analyse the following system from 8 
linear Diophantine equations 


l. aytan+a3 = S, 4. a,tast+a, = S, 7. aytastag = S, (5) 
2. agt+astadas = S, 5. a,+a5t+ag = S, 8. aastastaz = S. 
3. a7+ag+d 9 = S, 6. a3t+agtas = S, 


We note if symbols a, a2, ..., @9 are arranged as in the table 1, shown in figure, 
and their values are replaced by ones, which are taken from some total algebraic 
solution of the system (5), then table 1 will be transformed into the total algebraic 
formula of Magic squares of 3rd order. In other words, the discussed problem on 
solving the system (5) connects direct with the well-know ancient mathematical 
problem on constructing numerical examples of Magic squares of 3rd order. 


3.1 Requirements to a set of numbers, which is the solution of the system (5) 
Proposition 1. A set of nine numbers is a solution of the system (5) only in the case 


if one succeeds to represent these nine numbers in the form of such three arithmetic 
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progressions from 3 numbers whose differences are identical and the first terms of 
all three progressions are also forming an arithmetic progression. 


Proof. Using standard algebraic methods (for instance, Gauss method) we find 
that the total algebraic solution of the system (5) has the form 


a, = 2as—ag, a,=2a9+ag—2as; a3=3a5—dg—a9; a4 = 2a5— a6; (6) 
a7 =d9+dg—as; dg = 4as5—2a9— a6, 


where values of parameters as, ag and ag are chosen arbitrarily. Arranging solutions 
(4) in order, shown in the table 2 (see figure), we obtain the table 3. It is noteworthy 
that arithmetic progressions with the difference 2a5—a¢— ay place in the rows of 
the table 3, whereas ones, having the difference as —y, place in its columns. If one 
introduces three new parameters a, b and c by the equalities as=a+b+c, a=a+ 
2c and aj =a-+b into the table 3, then this table will acquire more elegant form, 
which it has in table 4, and so the fact of existing of the arithmetic progressions in it 
will receive more visual impression. Thus, the proof of Proposition 1 follows 
directly from the construction of tables 3 and/or 4 and it is appeared as a result of 
using the algebraic methods, mentioned in the points (a) and (c) of Sect. 1. 


"tap rar=tas [ass | 
Pag | Dasa | Aas — asa | 
(3) 


eS 
(4) 


Figure. Elucidating the general legitimacies existing between the elements of 
the solution (6). 


3.2 Elaboration of a universal algorithm for finding all numerical solutions of the 
system (5) from a given class of numbers 


Let it be required to find all numerical solutions of the system (5), which belong to 
the given class of numbers and has as=/f. For elaboration of a universal algorithm, 
solving this problem, we first write out all possible decompositions of the number 2f 
in the two summands of the following form 


2f = x(j)+ x2), (7) 


where j is the number of a decomposition; x(j) and x2(j) are two such numbers 
that x(j ) <x2(j ) and both ones belong to the given class of numbers. In a complete 
set of various decompositions of the kind, we fix only one, having, for instance, the 
number k. Determine for it the number d(k) 


d(k) = 2f- x(j). (8) 
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Proposition 2. The desirable numerical solution of the system (5), which con- 
tains as=f and numbers x,(k) and x2(k) of (7), can be found only in the case, if one 
succeeds to find, among the remaining numbers of the form x,(k), an arithmetic pro- 
gression from three numbers with the difference d(k). 


Proof. The truth of Proposition 2 follows from the construction of the tables 3 
and/or 4, shown in figure. 


It is evident also, to obtain a complete set of solutions of the system (5) from the 
given class of numbers, one should repeat the foregoing actions for all the differ- 
ences d(k). 


3.3 Deriving an analytical formula for counting the quantity of various solutions of 
the system (5) from natural numbers 


Proposition 3. Jf A(m) is the total number of various solutions of the system (5) 
from natural numbers and as=m_ then its value may be computed by the formula 


A(m) = 9[m/6}’ + (3 (m mod 6) — 8}[m/6] +2 — 2[{(m mod 6) +5}/6]+ (9) 
[(m mod 6)/5]. 


Proof. We first write out all possible decompositions of the number 2s in two 
distinct terms. 


2m = 1+(2m-1) (Gj=1, d(1) = m-1), (10) 
2m = 2+(2m—2) (j= 2, d(2) = m-2), 
2m = m—-2 +(m+2) (j=m-2, d(m-2) = 2), 
2m= m—-1 +(m+1) (j=m-1, d(m-1}1) = 1D). 


The problem on counting total number of various solutions of the system (5) with 
as =m is now reduced, in accordance with the universal algorithm of Sect. 3.2, to 
counting a total number of various arithmetic progressions consisting of three num- 
bers, which may be composed from the numbers 1, 2, ..., 7—2,m-— 1 and such that 
the differences in these progressions are respectively equal to d(m— 1), d(m-— 2), ..., 
d(1). 

To simplify this new problem we shall deduce a recurrence relation which will 
link the total numbers of various solutions having as=m and as=m-— 1. For this 
aim we decompose all the solutions with as;=m in two groups. The solutions, hav- 
ing number 1, will be attributed to the first group. A total number of such solutions 
will be denoted by Aj(m). All the remaining solutions we shall attribute to the 
second group. We decrease now each number by 1 in all solutions of the second 
group. After this operation, a lot of the second group solutions will represent by 
themselves a complete set of various solutions from natural numbers with as = m — 1. 
Thus, the following relation 


A(m) = Ay(m) + A(m=1) (11) 
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is valid or, in other words, if we know a value of A(m — 1) then for finding the value 
A(m) it will be sufficient to count the number of the solutions containing as = m and 
the number 1. This new combinatorial problem can be reformulated as the following 
one 


to find a total number of various arithmetic progressions from three numbers 
which can be composed from the sequence of numbers I, 2, ..., m— 2, m— 1 and 
such that the first number of these progressions is number I and the differences of 
the progressions are respectively equal to d(m—1), d(m—2), ..., d(1). 

It seems to be evident, that a total number of the desired progressions coincides 
with the maximal difference value of the progression Dmax for which one can still 
find an arithmetic progression of the required form from the set of numbers 1, 2, ..., 
m—2,m-—1. The value of D,,,, can be found from the correlation 1 + 2Da., =m-1, 
whence Drax = [(m — 2)/2], where square brackets denote the integer part. But in 
reality this value of Dmaox is not always coinciding with the value of A,(m): to elimi- 
nate this non-coincidence we must decrease the total number of arithmetic progres- 
sions by one if numbers 1 + d(k) or 1 + 2d(k) coincide with the number x,(k) of 
(7). 


Let us determine at which values of d(k) this coincidence occurs: 
1+ d(k,) = cic m— d(ki); 1+ 2d(k2) =k, =m — d(k2), (12) 


whence d(k,) = (m — 1)/2, and d(k2) = (m — 1)/3. If d(ki) = (m — 1)/2 the number 1 + 
2d(k,) > m — 1. Consequently, this case is never fulfilled. The coincidence occurs in 
the second case if m— 1 is multiple of 3. 

If we decompose all m-numbers in six classes so that the numbers of the form 6k 
will be attributed to the first class and those of the form 6k + 1 — to the second one 
and so on, where k = 1, 2, ..., then for all six classes of the m-numbers one can write 
out in the explicit form the values of Dax and Ai(?): 


m= 6k, Dmax = 3k-1, Aj(6k)=3k—1; (13) 
m=6k+1, = 3k-1, A,(6k+1)=3k-2; 
m= 6k +2, = 3k, A\(6k + 2) = 3k; 


m= 6k +4, = 3k+1, Aj,(6k+4)=3k; 
m=6k+5, Dmox = 3k +1, As(6k+5)=3k +1. 


Further we shall need the value of the difference AA(k, i) of the following form 
AA(k, i) = A(6(k +1) +i) — AGk+3), (14) 


D max 
Dax 

m= 6k +3, Dmax = 3k, Ai (6k + 3) = 3k; 
Dax 


where i = 0, 1, ..., 5. Using (14), (13) and (11) we may find an explicit expressions 
for AA(k, i). Let, for instance, i = 0. Then 


AA(k, 0) = A(6(k + 1)) — A(6K) = Ai(6(k + 1)) + Ai(6K + 5) + (15) 
+ Aj(6k + 4) + A,(6k +3)) + Ai(6k+ 2) + Ai(6k + 1) = 
= (3(k + 1)—1} + Gk + 1) + 3k + 3k + 3k + (3k-2) = 18k +1. 
Remaining values of AA(k, i) for i = 0, 1, ..., 5 can be found analogously: 
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AA(k, i) = 18k + 1 + 3i. (16) 


It is evident, since the AA(k, i) is a linear function from k, values A(6k +i) may 
be obtained from the second degree polynomial 


A(6k +i) = boi) + by) k + boli), (17) 


where b2(i) =9, bi(i) = 31-8, boi) = 2—2[(i + 5)/6] + [W/5]; square brackets mean 
the integer part. Taking into account that i=(mmod6), k=[m/6] and 6k+i =m 
we may obtain (9) from (17). 

It should be noted that using regression analysis methods one may appreciably 
simplify *° the expression (9): 


A(m) = g{(3m? — 16m + 18.5)/12}, . (18) 


where the notation g{a} means the nearest integer to a. 


4 Algebraic approach to deriving analytical formulae available for 
constructing classical Magic squares of the m-th order 


We remind that in the general case’? Magic squares represent by themselves 
numerical or analytical square tables, whose elements satisfy a set of definite basic 
and additional relations. The basic relations therewith assign some constant property 
for the elements located in the rows, columns and two main diagonals of a square 
table, and additional relations, assign additional characteristics for some other sets of 
its elements. In particular, when the constant property is a significance of sum of 
various elements in rows, columns or main diagonals of the square, then this square 
is an Additive one. If an Additive square is composed of successive natural numbers 
from 1 to n’, then it is a Classical one. 

It is evident*’*"* that, from the point of view of mathematics, the analytical 
solution of the problem on constructing Classical squares of the n-th order consists 
of determining a form of fand g functions, which permit to compute the position for 
any natural number N from 1 to n’ in cells of these squares: x = f\N, n) and y = 9(N, 
n). 

In this section we 

1) adduce two types of analytical functions, by which one may construct Classi- 
cal squares of odd orders; 

2) reveal a connection between these analytical functions and Latin squares; 

3) give an algebraic generalisation of the notion “Latin square”; 

4) derive analytical formulae available for constructing Classical squares of both 
odd and even orders. 


4.1 Classical approach to deriving analytical formulae available for constructing 
Magic squares of odd order from natural numbers 


For any linear algorithmic methods of constructing Classical squares of odd orders, 
the functions f and g have the following forms ***'': 
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AN, n) = a(N-1) + b[ (N= Din] + c (mod n), (19) 
(WN, n) = ax(N-1) + b[ (N- lin] + co (mod n), 


where square brackets mean the integer part; a sign “S” is the modulo n equality; N 
is any natural number from 1 to 7; a), bj, c, and a, by, cz are such integral coeffi- 
cients, that the numbers aj, a2, b;, bo; a,b. —a2b,; a2—a1, )— bi, a2 +a, and + 
b, are mutually disjoint with n. 

There is no difficulty in counting that in formula (19) the coefficients {a), bi, ci, 
do, ba, C2) are equal > 10, 11 

— { 1, 1, —[n/2]; 1, -1, [n/2] } for Terrace algorithmic method of constructing 
Classical squares; 

— {1,-1, n/2]; 1,-2, n — 1} for Siamese method; 

— {1,-1, [n/2]; 2, 2,0 } for Knight method; 

— {(3-a)q+ (a+ 1)/2, 3 -—a)q+ (a—1)/2,0; (3 -a)q + (a - 19/2, (3 - a)g 
+ (a+ 1)/2, 0 } (where q = [(n + 1)/6], a =n —6q) for the classical square of the n- 
th order, which, if n is an odd number, non-divisible by three, can be formed also 
from a pair of orthogonal Latin squares, constructed by the pair of comparisons x + 
2y (mod n) and 2x + y (mod n); and so on. 

It should be noted, that the above conditions for coefficients of functions f and g 
become contradictory for even n. For example, by the conditions, the coefficients 
@\, G2, b, and b, of the functions f and g of (19) should be mutually disjoint with n, 
and consequently, if n is even, they must be odd. The same requirement must be the 
true for the number d = a,b — apb;. But if a, a2, b; and by are odd, the number d, 
which is the difference of the two odd numbers, will be an even number. 

Thus, an essential fault of linear formulae of (19) is the impossibility of using 
them for constructing Classical squares of even orders. 


4.2 Revealing a connection between Latin squares and analytical formulae of (19) 


Proposition 4. [f a Classical square of the n-th order is constructed by formulae 
(19), then it may be constructed also by the formula 

Nix, y) = n p(x, y) + r(x, y) + 1, (20) 
where p(x, y) = %x+ Biyt+ oO; and r(x,y) = Gxt Pry + O2. 

Proof. The equivalence of formulae (19) and (20) appears from their linearity 
and the fact, that (20) are inverse formulae to (19). In particular, if values of 


coefficients (a), b1, C1, @2, b2, c2} of formulae (19) are known then values of {04, Bi, 
61, Q, Bo, 02} of formulae (2) may be computed from following linear equations " 


mo, = — a2; mB, = ay; mo, = anc, — arco; Cy 
m2 = bo; mB. = — bi; mo2 = bic2 — boc; 
m = a\b2 — and, 


and, reciprocally, at the reverse task, the values of {a), bi, ci, @2, b2, cz} may be 
computed from equations: 
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Ha;= ~Bi; pbi= Bo; pc: = Bio2 — foou; (22) 
Ha2= 0); Hb, = — 02; Ue2 = O20) — 0102; 
w= ai Be a 2B. 
For instance, values of {01, B:, 01, G2, Bo, 62} of formulae (20) are equal 
— {(n+1)/2, (n-1)/2, (n—1)/2; (n+1)/2, (n+1)/2, 0 } for Terrace method; 
— {1,n-1, (n—1)/2; 2, n—1, (n—1)/2 \} for Siamese method; 


— {[n/2], a, n-a ; (n+1)/2, a, a } for Knight method, where a = dc + (n— d) (1- 
c), d=[(n+1)/4], c = [(n mod 4)/2]. 


We remind* that, a quadratic table nxn in size is Latin square of n-th order if 
only n elements of this table are different and each of these n elements occurs only 
one time in each row and column of the table. The two Latin squares P and R of the 
same order n are called orthogonal if all the pairs, formed by their elements p;; and 
r;j (i is the number of a row; j is the number of a column) are different. 


Proposition 5. [f elements of a Latin square of the n-th order are numbers 0, 1, 
... n= I, then, for constructing such Latin square, one may use a linear comparison 


L(x, y)= ax+ By+ o, (23) 


where a@ and B are integer numbers, which are to be mutually disjoint with n; o is 
any integer number. 


Proof. Let the numbers L(x, y) of (23) are located in each cells of a quadratic 
table nxn in size. We consider nm numbers, which are located in row yo of this table. 
Since the discussed numbers are obtained from the linear comparison (23) at x = 0, 
1, ..., 2 — 1, to show that all they are different, we should demonstrate that they 
belong to different modulo n classes. Let x1 > x2 and ax+ By + GO = ax 
+ Byo +o. Since B yo+ o is a constant, in accordance with the properties of com- 
parisons **, we obtain the new equality ax; = x2. Hence, since « is mutually dis- 
joint with n, x, = x). But this equality contradicts our assumption. Thus, each of 
numbers 0, 1, ..., 2 — 1 occurs only one time in each row and column of the dis- 
cussed table and so this table is Latin square of n-th order. 


Proposition 6. Every Classical square of the odd order, decomposed on two 
orthogonal Latin squares, may be constructed by the formulae (19) and otherwise. 


Proof. The truth of Proposition 6 follows directly from Propositions 4 and 5 and 
conditions for coefficients of functions f and g of (19). 


4.3, Deriving analytical formulae available for constructing Classical squares of 
both odd and even orders 


The way 1. Let us give an algebraic generalisation of the notion “Latin square”: 


a quadratic table nom in size is the generalised Latin square of n-th order if only 
n elements of this table are different and each of these n elements occurs only n 
times in this table. 
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Proposition 7. Every Classical square of a order n may be decomposed on two 
orthogonal generalised Latin squares P and R of the order n. 


Proof. To prove Proposition 7, it is sufficient to note that 
a) any integer number NV from | to n? may be presented in the form 


N=np+r+l, (24) 


where p and r can take values only 0,1, ..., —-1; 

b) each of values 0,1, ...,—1 of parameters p andr occurs n times precisely 
in the decomposition (24) of numbers N. 

Thus, to construct two orthogonal generalised Latin squares P and R from a 
Classical square of a order n, one should replace in the Classical square all numbers 
N by respectively (N-—1) modn and [(N-1)/n]. 


Proposition 8. Every Classical square of order n may be constructed by the 


formula (20), in which functions p(x, y) and r(x, y) may belong, in general case, to 
both linear and non-linear classes of ones. 


Proof. The truth of Proposition 8 follows directly from Propositions 7 and mate- 
rials of Sect. 4.1. 


We note, in particular, one may construct Classical squares of even-even orders n 
(n = 4k; k = 1,2, ...) by the analytical formula (20), in which functions p(x, y) and 
r(x, y) have the following forms ieee 


Px, y) = cx + (1-c) (n-x-1) and r(x, y)=(l-c) y+c(n-y-1), (25) 
where c = {[ (x+1)/2 ] + [( y+1)/2]} mod 2; or 
p(x, y) =cd—x—1+(1—c)(n—-d) and r(x, y)=by+ (1-5) (n-y—-)), (26) 


where c = (x + y+a)mod2; d=(l-a)y +a(n-—y-1); b= {[(x+3)/2] + [ y/2] + 
a} mod 2; a=[2y/n]; and so on. 


The way 2. It is evident, we may consider Classical squares not only as the sum 
of two orthogonal generalised Latin squares (see the way 1) but, for instance, as 
quadratic tables whose rows contain certain numerical sequences. Let us look into 
the problem on finding universal analytical formulae for constructing Classical 
squares from this new point of view. 


Proposition 9. If a Classical square of the n-th order is constructed by formulae 
(19), then it may be constructed also by the formula 


Na y= at+b-Ac, (27) 
where a,bandc are any integer numbers; A is 0 or 1, the sign “=” is the modulo 
n’ equality. . 


Proof. Let a Classical square of an odd order is constructed by formulae (19). It 
follows from Proposition 4 that this square may be constructed also by formulae 
(20). We deduct x-th element of first row from every x-th element of all y-th rows of 
the Classical square. It is evident that the number 
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{ n{(a, x + Bi y + 6;) modn} — n{(a, x +01) mod n} + (28) 
(Q2.x + Boy +62) modn — (2.x +62) modn ) mod n’. 


will be located in the cell (x, y) of the reformed Classical square {see (20)}. Using 
the equality (dn) mod n? =n (d mod n) we present (28) as the sum of two summands 


n{(Bi y)} mod n + {(O2 x +Bo y +62) mod n — (2 x + G2) mod n}. (29) 


The second summand of (29) may have only two values: (B2y) modn or n-— (fp y) 
mod n. Thus, we obtain that numbers of any y-th row of the reformed Classical 
square may have only two values. By using the mentioned method of constructing 
formula (27), we find that parameters of this formula a, b, c and are connected 
with parameters of the formula (20) by correlations 


a= n{(x+01) modn} + (Q)x+62) modn, (30) 
b =n((Bi y) mod n} + (By) modn, c=n, 
XA ={[1-— sign{(Q2x + fo y + (64) modn— (a) x + (64) modn } 1/2, 

where sign(x) =|x Ix if x#0 and sign(0) = 0. 


It should be noted, if we get off the sign “=” in the formula (27) and translate 
correlations (30) into language of numerical sequences {see the point (c) of Sect. 1}, 
we obtain that, for algorithmic methods which mentioned in Sect. 4.1 and 4.2, the 
parameters of formula (23) are determined by correlations 


a =—(-1)' n(n— 19/4 + k(n + 1/2 +(n(n-3) +2}/4, (31) 
b=n-1-y, c=n, A=[sign{(h—-y) +2}/2], 
h=[z/2]-1+(n+1)(z2-[z2]), 

0,(z)=y+z+2—n[sign{ (y+z—n+1)+2 }/2], ki = 0,0) 


where the numerical sequence {a}, if its values are computed atk = 0,1, ...,.n—-1, 
coincides with the numerical sequence, located in the first row of the Classical 
square; 6,(t) is a permutation operator of numbers 0, 1, .. ., 1 — 1 and for 


— Terrace method k=k,, z=k,—1; 

— Siamese method z= (n—2(kj-—1)} modn, k= z+; 

— Knight method ky = 6,(ki- 1), z= —(-1)" n/4 + kp/2 + (nt4)/4, k= 241. 

It is evident, using the formula (27) with parameters (31), one can have no diffi- 
culty in discovering “genetic connections” between different Classical squares 
and constructing methods and in generating a set of new methods. For instance, if n 
is an odd number, non-divisible by three, the new algorithmic methods for 
constructing Classical squares of odd orders appear when k, = oy (x) or k, = Gy (x) 
in (31), or the form of 6, and/or the numerical sequence {a;} is changed. 

It remains for us to add that parameters of the formula (27) are determined by 
correlations 


a=nk, b=w, c=n-2w-l, X=[((k+2) mod 4)/2], (32) 
Oz)=1+{z—- h(2((z+h) mod2)-1)} modn, 
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h=y+c[yly/2]], k=o.(z) 


for formulae (25) and (26), where oy (z) = 1 + {z—A (2((z + A) mod 2) -1)} mod n; h 
=y+c[y/[ y/2]]; k, =0,(z) and for 
— the formula (26) z=x, w=y; k=k,; 


— the formula (25) 4, = [((y + 1) mod 4 )/2]; k= ky; A, = [((x+1+2A,) mod 
4/2]; z =(4;+n—h(1—-2(x, mod 2))) modn; x, =Agxt+(l-A) (n-x-1l); 
w=Ayy+1—A.)(n-y-l). 
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